Charge-induced instability and macroscopic quantum nucieation phenomena at the 

crystal 4 He facet 
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An existence of the charge- induced instability is well known for the 4 He crystal surface in the rough 
state. Much less is known about charge-induced instability at the 4 He crystal surface in the smooth 
well-faceted state below the roughening transition temperature. To meet the lack, we examine here 
the latter case. As long as the electric field normal to the crystal facet is below the critical value 
same as for the rough surface, the crystal faceting remains absolutely stable. Above the critical field, 
unlike absolutely unstable state of the rough surface, the crystal facet crosses over to the metastable 
state separated from new crushed state with a potential barrier proportional to the square of the 
linear facet step energy. The onset and development of the instability at the charged crystal facet 
has much in common with the nucieation kinetics of first-order phase transitions. Depending on 
the temperature, the electric breaking strength is determined either by thermal activation at high 
temperatures or by quantum tunneling at sufficiently low temperatures. 

PACS numbers: 67.80.-s, 64.60.Q- 



5-H 



C3 



-a 
c 

o 
o 



> 

o 
o 

(N 



X 



I. INTRODUCTION 

It is well known that a charged interface between two 
fluids can develop an electrohydrodynamic instability at 
sufficiently high density of charges. Such charge-induced 
instability results from the competition between the elec- 
tric forces of like charges and forces of surface tension 
and gravity. Naturally, the liquid phases of helium have 
become one of physical systems for the theoretical and 
experimental studies of electrohydrodynamic instabili- 
ties in particular, as softening of the gravitational- 
capillary wave spectrurn[2|-y] , charge- induced deforma- 
tions of the interface [5|-l7| , formation of regular array of 
dimples 0-0], individual multielectron dimples 0, and 
hexagonal reconstructuring @. 

It is of particular interest to compare the onset and de- 
velopment of the electrohydrodynamic instability at the 
liquid-solid 4 He interface with that at the interface be- 
tween two fluids. The first theoretical and experimental 
studies have shown that a charged-induced instability at 
the superfluid-solid 4 He looks roughly like the electro- 
hydrodynamic instability at the free liquid 4 He surface 
[10l4l3l | . For the high temperature region where the crys- 
tal surface is in the rough state, such similar behavior 
is expectable since the superfluid-solid interface in the 
rough state has a very high mobility and interface exci- 
tations represent weakly damping crystallization waves 
whose dispersion [l4|, HH is quite similar to that of usual 
gravitational-capillary waves at the free liquid surface. 

To date, no systematic study has been made on the on- 
set and development of electrohydrodynamic instability 
at the well-faceted and atomically smooth crystal sur- 
faces which may have an infinitely large stiffness and 
excitation spectrum differing from the usual crystalliza- 
tion wave spectrum. The most striking distinction of the 
smooth faceted crystal surface from the rough one is the 
existence of nonanalytic cusplike behavior in the angle 
dependence for the surface tension, e.g., [l6l - [l8l |. The ori- 



gin of the singularity is directly connected with nonzero 
magnitude of the facet step energy below the roughening 
transition temperature of about 1.2 K. 

In present work we attempt the electrohydrodynamic 
instability at the smooth faceted surface of a 4 He crystal 
in contact with its liquid phase. As we will see below, the 
close similarity between the rough and smooth states of 
the crystal surfaces extends until the charge density is be- 
low the critical one and the state and shape of the surface 
are stable. As the charge density increases, the develop- 
ment of the instabilities becomes different in kind. Un- 
like the rough crystal surface, the faceted surface crosses 
over to a metastable state and the further development of 
the instability is accompanied by overcoming some nucie- 
ation barrier. The barrier height is proportional to the 
square of the facet step energy and drastically reduces 
as the charge density increases. At the sufficiently low 
temperatures the thermal activation mechanism of over- 
coming the barrier is replaced with the quantum tun- 
neling through the nucieation barrier. On the whole, 
the charge-induced reconstructuring of the faceted crys- 
tal surface resembles much first-order phase transitions 
and macroscopic quantum nucieation phenomena in the 
helium systems |19l - l23j . 

For simplicity, we keep in mind the basal plane of 
hexagonal 4 He crystal as an example of the crystal facet 
and neglect any anisotropy in the plane. We also suppose 
that the temperature is below the roughening transition 
temperature and the crystal surface is well-defined and 
faceted. 



II. HAMILTONIAN. THE ONSET OF 
INSTABILITY AT THE CRYSTAL FACET 

Let us assume that the crystal surface is parallel to the 
xy plane, with vertical position at z — 0. In order to con- 
sider the stability of the surface, we proceed as follows. 
First, we call C = C( r ) the displacement of the surface 
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from its horizontal position z — with r — (x, y) as a 
two-dimensional radius- vector. In addition to the surface 
tension force and the force of gravity due to difference in 
the densities between the solid and liquid states Ap, one 
should involve also the interaction of the charges with 
electric field E and the direct Coulomb interaction be- 
tween the charges. Then the total energy U of a charged 
surface can be written as 



U 



d 2 r (a(u) y/l + (V() 2 + Apg^- + eEn(r)( 
, en(r)en(r') 



d 2 rd z 



\l-V 



(1) 



Here I = {r,Q stands for the three-dimensional coor- 
dinate of a point at the surface, n(r) is the density of 
electrons with charge e, g is the acceleration of gravity, 
and a(v) is the energy of a unit surface area or surface 
tension. 

Unlike the fluid-fluid interface, the surface tension co- 
efficient a(v) for the crystal facet depends essentially on 
the direction of the normal v to the interface against crys- 
tallographic axes. In our simplest case this is a function 
of angle t? alone between the normal and the crystallo- 
graphic [0001] or c-axis of the crystal hep structure with 
the geometric relation | tan$ | = | V£ |. 

For the crystal facet tilted by small angle d from the 
basal plane, the expansion of surface tension a(fi) usually 
written, e.g. |16l - [l8| . as 

a(i?) = (a + aitan | i? | +...)cos0, | tantf |=| VC |, 
can be represented for the small angles by a series 
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We intentionally do not write the next terms of expan- 
sion, e.g., cubic one due to step-step interaction, since 
we assume to study only small bending of the crystal 
surface. The angular behavior has a nonanalytic cusp- 
like behavior at $ = due to a\ = a\{T) representing 
a ratio of the linear facet step energy j3 to the crystallo- 
graphic interplane spacing. Below the roughening tran- 
sition temperature for the basal plane Tr ~1.2 K the 
facet step energy j3 = P(T) is positive and vanishes for 
temperatures T > Tr. 

To determine the equilibrium shape of the surface £(r) 
and equilibrium charge distribution n(r), we must mini- 
mize the energy functional (TTJ) against C and n at a given 
total surface charge Q. This condition can readily be 
taken into account by augmenting the energy functional 
with the Lagrange factor A in the form 

-A J en(r)d 2 r. 

In addition, treating the energy functional, we naturally 
imply one more obvious condition n(r) 0. 



In the general form the minimization of the energy 
functional is a practically unsolvable problem because 
its expression ([!]) contains not only quadratic terms in C 
and n. Thus, we start first from analyzing small surface 
bending £(r) and small gradients |V£| -C 1. The latter 
implies implicitly that |C| -C r and we can put approx- 
imately I = r in the denominator of the Coulomb term 
in Eq. ([T]). Next, we expand the surface tension term in 
| VC|j retain the linear and quadratic terms alone, and ar- 
rive at the following expression for the total excess energy 
U associated with nonzero surface bending ((r) 
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Here we have labeled 
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a = a(0) +a"(0) 

as a surface stiffness. 

The spatial scale of surface distortion is usually deter- 
mined by the capillary length Ao = (ao/Apg) 1 / 2 w 1 mm 
if one takes ao ~ 0.2 erg/cm 2 and Ap s» 0.018 g/cm 3 for 
4 He [l6|. The electric field E and charge surface density 
en have the same dimensionality and their typical scale 
is (aoAp g) 1 / 4 w 400 V/cm. Accordingly, the typical 
electron density equals {a^Ap g) 1 / 4 / e ~ 2.8 x 10 9 cm~ 2 . 
The number of electrons ira \Apg)-3/ A within the cir- 
cle of radius Ao runs to 10 8 . And lastly, unit of energy 
corresponds to a^/(Apg) sa 2.2 x 10 -3 erg. 

So, if we measure Q and r in units of capillary length Ao, 
electric field and charge density in units of {a^Ap g) 1 ^ 4 , 
and energy in units of a^/ {Ap g), the total excess energy 
U can be expressed in terms of dimcnsionless units as 



U = 



d 2 , 



a 



(VC) 2 



L + E n(r)<: 



d 2 rd 2 r> n(r)n(r>) 
\r — r \ 



(2) 



As for the step energy a%, we assume that its low tem- 
perature value [lH is approximately ct\ f=s 0.014 erg/cm 2 . 
This value amounts to one-tenth of surface stiffness ao 
and in the following we can keep inequality ai/«o *C 1 
in mind. Moreover, this small parameter justifies all ap- 
proximations that will be made further. 

The uniform state of the surface holds for the elec- 
tric field values as long as the contribution to the excess 
energy ([2]) due to variations C( r ) from C = and n(r) 
from homogeneous value h is a positive-definite quantity. 
Using the following equation for Lagrange factor 



and putting 8n(r) = n(r) — n, we find for the variation 
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of the excess energy 

ST' = ld 2 r(^\V(\ + 
1 

+ 2 



(vc) 2 , e 



E5n(r)( 



2 2 , 5n{r)Sn(r') 
a r a r — ; — . 
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To analyze it, we use the Fourier representation 



C(r) = Y,(ke ikr and Sn(r) = ^<5n fe e 

k k 

and rewrite the energy variation as 



ikr 



su 
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l)CfcCfc + E(Sn k Q + Sn* k ( k ) 
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Minimizing SU over Jn^ yields the optimum relation 

k 

Sn k = -— Cfc (4) 

Z7T 

and the corresponding optimum value of the energy 

5U = [ ^|VC|d 2 r+i^(fc 2 + l-^)|C fc | 2 . (5) 
i do 2 fc 2?r 

The second term is always positive provided the inequal- 
ity E 2 < 2ir(k + 1/k) is satisfied for all wave vectors k. 
The minimum of the right-hand side of the inequality oc- 
curs at k = fe c = 1 and corresponds to the critical field 
E c = Thus, the crystal facet is absolutely stable at 

E < E c . 

At E > E c the stability is lost and the distortions of 
the homogeneous state should appear. In this regard the 
situation resembles the loss of stability for the rough state 
of the crystal surface. However, the development of the 
stability and the transition to unhomogeneous state differ 
drastically. In fact, due to positive a± > term linear 
in |VC| the evolution of the crystal facet perturbations 
should inevitably be accompanied with overcoming some 
potential barrier, the barrier height being dependent on 
the field strength E. The more the field strength, the 
less the potential barrier height. 

To proceed, let us return to the coordinate representa- 
tion of Eq. fl5J) 



SU 

1 

~2 



ao 2 2 

, E 2 (V r C(r)V r ,C(r')) 

a r a r 



(6) 



and give a qualitative description of the matter. For this 
purpose, we employ a variational principle and dimen- 
sional analysis of the functional ©. Let us represent the 
surface distortion £(r) with the aid of the trial function 
f(x) in the axially symmetrical form as 



C(r) = (f(r/R), 



(7) 



where £ is a typical magnitude of distortion and R is its 
typical size. Then we have 

and the dimensionless factors are given by 



o= / \f'(r)\2nrdr, b= I f' 2 (r)2Trrdr, 
Jo Jo 

c = J™ f 2 (r)2Trrdr, d = J d 2 r d 2 r' T ' — f '^ f '^ 



OO / pOO 



dk 



2 



dr2nrf'(r)J 1 (kr) 



where J\{x) is the Bessel function of the first kind. 

As the electric field strength exceeds the value E = 
(87r 2 v / 6c/d) 1 / 2 , there appears a region of £ and R with 
the negative values of SU separated always from SU — 
at C = with the intermediate positive SU values. 
Rewriting the excess energy 6U(£, R) as 



SU=- 
c(E 2 



lf8n 2 ai \ 2 a 2 b/d? 



Of) 



E± - E<k 



c E 2 



bE% 



ai 



a/c 



a E*/E% ~ 1 



E 2 



in 2 Vbc 



one can readily see that the state of the crystal facet 
changes from the stable to metastable state at E > E 
and the point 
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and Rq 
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E 2 



becomes a saddle point of the potential relief. The po- 
tential barrier height equal to 



•2 '2 7714 

a a x E 

UO — n 



2c a 2 G E 4 - E 4 

must be overcome to break the flat faceting of a crystal 
surface. 

Unfortunately, we cannot find the exact function f(r) 
and, correspondingly, values of factors a, b, c and d which 
optimize the functional ©• However, it is clear that the 
potential barrier height should be infinitely large at E = 
E c and thus Eq = E c for the exact solution. This entails 
the obvious relation d — 2n(bc) 1 / 2 between coefficients 
for the exact solution. To estimate them, we use a trial 
function f(x) = exp(— x 2 ). The direct calculation results 



and 
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2nVbc 
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0.94 



(8) 
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in place of unity. Hence we may expect an accuracy of 
our estimate within about 10%. 

Let us compare the height Uq of the potential barrier 
at the saddle point with the roughening transition tem- 
perature Tr about 1.2 K. In the dimensional units we 
have 



U 



2c Apg E 4 - £4 



1.4 x 10 11 



Et 



E A - Et 



(in K). 



One may be surprised with the huge barrier height so 
that, unlike the rough crystal surface absolutely unstable 
at E ^ E c , tens of E c should keep a crystal facet prac- 
tically stable for an experimentally available time. Pro- 
vided we expect a reasonable observation time of destruc- 
ting the faceted state due to thermal activation mecha- 
nism, we should provide a ratio Uq/T of about a few tens 
[2l|,|22j This means that the electric field E should exceed 
the critical one E c by a factor of about 300. The same 
factor certainly refers to the surface density of charges. 

In the dimensional CGSE units the bending deflection 
Co and the typical size of inhomogeneity Rq are given by 



ia i , a E 
K ° ]=i7: Wc ai E^ 



i?0 



Et 

b 47r«o 
c E^^ 
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In the weak fields of few critical values the critical pa- 
rameters i?o and | Co | prove to be of macroscopic sizes 
in accordance with macroscopically large height of the 
potential barrier. 

For E = 300E C , we find approximately i?o ~ 16 nm 
and | Co I ~ 2 nm. On the whole, the electric field should 
be very large compared with the critical value E c in order 
to reduce significantly the nucleation barrier for the ef- 
fective production of a few circular crystal terraces tilted 
with the angle about arctan(cti/ao) ^4°. In this sense 
the critical fluctuation represents a region of the crystal 
surface in the rough state. 

From the physical point of view the angle of slope 
arctan(ai/ao) ~ 4° is determined by the competition of 
two contributions into the total surface energy. One orig- 
inates from the regular surface term «oC 2 and the second 
does from irregular step tension term ai|C|-R- Provided 
a ( 2 3> ai|C|-R, the latter contribution becomes negligi- 
ble and thus the interface properties resemble those in the 
rough surface state. On the contrary, if a C 2 ai|C|-R, 
the dominant term linear in |C| is responsible for the ori- 
gin of a potential barrier since the other terms quadratic 
in C are yet insignificant. 

Note that the small gradient approximation we use is 
satisfied since \\7(\ ~ |Co|/-Ro ~ oli/olq <C 1 with the 
exception of narrow region E ~ E c . The latter remark 
refers also to justifying small density variations 8n <C n 
valid to the extent of smallness |Co|/-Ro- 



III. LAGRANGIAN. THE QUANTUM 
BREAKING OF THE CRYSTAL FACET 

The destruction of the faceted crystal surface is accom- 
panied by overcoming some potential barrier depending 
on the charge surface density. There are two basic mech- 
anisms to overcome the potential barrier. One is the 
thermal activation efficient at high temperatures and the 
second is the quantum tunneling through a potential bar- 
rier dominant at sufficiently low temperatures. In order 
to treat the quantum tunneling, it is necessary to involve 
the interface dynamics, in particular, to determine the 
kinetic energy of the charged interface in addition to the 
potential energy U. 

As a first step, we employ the so-called metallic ap- 
proximation. In this approximation it is assumed that 
the mobility of electrons along the supcrfluid-crystal He 4 
interface is very high and the charged helium interface 
represents an equipotential surface so that the electric 
field is always normal to the interface as for a well- 
conducting metal. A necessary condition for such ap- 
proximation assumes at least that the plasma oscillation 
frequency of a two-dimensional layer of electrons with 
effective mass m e 



(2Trne 2 k/m e ) 1/2 



(9) 



is much larger than the typical frequency u> of the 
gravitational-capillary or melting-crystallization waves at 
the same wave vector k. So, within our first approxi- 
mation we believe that the charge density distribution 
n(t, r) has sufficient time to accommodate to the surface 
distortion ((t, r) and is determined by the electrostatic 
relations in accordance with the profile £(t,r). 

Neglecting possible energy dissipation, we describe the 
charged interface dynamics using the following action 



S 



dtL[((t, r),C(t, r),n(t,r)] 



(10) 



with the Lagrangian L equal to the difference between 
the kinetic energy functional and the potential energy 
functional U introduced by Eq. |T]) 



L = ^5 / / d~nrr 



i„ A 2j C(*, r)C(t, V) 



27r|r - 



-U[((t, r),n(t, r)}. 



Here we ignore the compressibility of the both liquid and 
solid phases. Because of low temperature consideration 
we will also neglect the normal component density in the 
supcrfluid phase or, equivalently, difference between the 
superfluid density p s and the density of the liquid phase 
p. Then the effective interface density p c g is given by 

Pes = (p' - pf/p ~ 1.9mg/cm 

and depends on the difference Ap = p' — p between the 
solid density p' and the liquid density p. For our pur- 
poses, the exact magnitude of the effective density is 
inessential. 
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Next, for convenience, let us introduce units of time 
equal to (/Ocff-^o/c^o) 1 ^ 2 ~ 3.1 ms and measure the action 
in units of (a peffAo) 1//2 ~ 0.62 x 10~ 5 erg-s. Using the 
speculations and arguments bringing us to Eq. Q and 
then to Eq. (|6]), we arrive at examining the following 
effective action 



S = J dtL eS [((t, r),C(t, r)} 



with the dimensionless Lagrangian 



L v tf — — 



C(t, r)C(t, r') 



d z rd z 



d 2 r[ — |VC| 



2n\r 
(VC) 2 



( 2 



, 2 , 2 , £ 2 (VrCft,r)Vr<(*,r / )) 
r (2tt) 2 |r-r'| 



Within an exponential accuracy the quantum decay 
rate of the metastable state is proportional to 

r oc exp(— Se/K), 

where Se is the effective Euclidean action calculated at 
the optimum escape path. This path starts at the en- 
trance point under the potential barrier and ends at the 
point at which the optimum fluctuation escapes from the 
barrier [2ll l22j |. In other words, quantum fluctuation 
penetrates through the potential barrier along the path 
of least resistance. Before calculating the quantum rate 
at which the crystal facet breaks up, we must go over to 
the effective Euclidean action defined in imaginary time 
t — > it. We refer to books [24], [25| for details. 

As a result, we should analyze the following functional 
defined within the time interval [-K/2T, K/2T] 



1 



dtL E [((t, r),C(t, r)], 



p rdV C(*,r)C(t ,r') 
2n\r - r'\ 

2 a2 



2., 

a 



(VC) 2 C 
2 2 



d 2 rd 2 r 



, E 2 (V T .C(*,r)Vr<(t,r / )) 
(2tt) 2 \r-r'\ 



Again, the exact determination of extrema for the action 
Se is a rather complicated problem. We here consider 
only the case of zero temperature when the limits of in- 
tegration over imaginary time are infinite. As before, it 
is convenient to take an advantage of the dimensional 
analysis and variational principle. We will express the 
surface distortion ((t,r) in the terms of function f(y,x) 
with the scaled arguments as 

t(t,r) = {f(t/T,r/R). 



Next, we calculate the action Se at zero temperature 



S E (C, t, R) = F^- + r ( ^A\C\R + B C - 
It \ao 2 



C 2 R 2 



-C 



D 



E 2 C 2 R 



2 4tt 2 2 

The numerical factors are given by the integrals 
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A = 
B = 
C = 



D 



dt / \f'(t,r)\2%rdr, 

oo JO 

dt / f 2 (t,r)2irrdr, 

oo JO 

/•OO 

dt / f 2 (t,r)2irrdr, 
Jo 

,r-r> f(t,r)f>(t,r>) 



dt / d 2 rd 2 r 



dt / dk 

J 

o 



dr2nrf'(t,r)J 1 {kr) 



F = 



dt / d 2 rd" 



r , f(t,r)f(t,r') 
2ir\r — r'\ 



dt 



dk 
2^ 



dr 2nrf(t, r)Jo(kr) 



where Jo(x) and J\{x) are the Bessel function of the first 
kind. 

From the condition of vanishing derivatives in £, R and 
r for Se we find the following parameters of the quantum 
critical fluctuation 



iCal 



Ra 



A a x 

7 [B_E% 
~ 4 V C E 2 



1 



E 2 E 2 



1M + 

16 E 4 4 J E 4 - E, 



T_E% 
16 E 4 
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1/4 



E 3 



1 - 



16 -E 4 4/ 



2 v / 2\ C 3 / E^E 4 - E i 



1 - 



16 E 4 - 
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Here E = (8tt 2 VBC / D) 1 ' ' 2 which should coincide for 
the exact solution with the critical field value, i.e., Eq = 
E c = v4~7r. Then we calculate the corresponding value 
of action S q according to 

Sq = — \Cq\RqTq 

2 ao 

at the critical point R q ,r q ) representing a saddle 
point of the functional S E (fTTj) . Finally, we obtain 



^ 2 )^ ai 

q 16V2 CV 4 a 2 ' 



49A 2 {BF 2 



1 _ JL B o 

x 16 W 



i) 
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x 16 W 



En 



E \E 4 -E 4 
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Like the potential barrier height, the action S q becomes 
infinite at the same critical field E = E . 

To estimate the numerical coefficients F, A, B, C, 
and D, we choose a physically expedient trial function 
f(t,r) — exp[— {t 2 + r 2 )]. The straightforward calcula- 
tion gives 

7T 3 / 2 /tt\ 3/2 7T 3 7T 2 

2V2' V 2 / 2 ' 4 

with the same ratio E c /E as in ©. 

Let us compare the action S q with the Planck con- 
stant h. Introducing a facet capillary length Ai = 
(qi/ Ap g) 1 ! 2 , we have in the dimensional units 



Sq = 49A 2 (BF 2 ) 1 / 4 f^d^ PcS Xl 
h ~ 16V2 C 7 / 4 



h 2 



i L £5 

x 16 W 



lf 2 E c 



E\E 4 - E 4 



1/2 



3/2 



16 IF 



4 x 1CH 



_7_ -B* 
16 IF 
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f 
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2 E \E 4 - E 4 
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As is seen, even for the electric fields which are dozens of 
times larger than the critical one E c , the ratio S q /h has a 
giant magnitude so that the crystal surface will remain in 
the well-defined faceted state for the practically infinite 
time. At E ^> E c we have an estimate 

S q /hn 2.3 x 10 21 {E C /E) 7 . 

For strong E 3> E c fields, in the dimensional CGSE units 
the bending deflection £o and the typical size of inhomo- 
geneity Rq are given by 

1*1 . 7A a\ 

ICgl = 47T- 



Oil 



7 B 47rao 
9 _ 8 V C £ 2 



15^2. 



For E = 300i? c , we find approximately R q ~ 13 nm and 
~ 5 nm. The estimate of the tunneling time in the 
strong E ^> E c fields yields 

^ 2 {E C /Ef (in seconds). 
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Again, the small gradient approximation is fulfilled 
since |V£| ~ \C, q \/R q ~ ctx/ao *C 1. Let us compare 
the plasmon frequency VL p with the inverse time of tun- 
neling t~ 4 in order to justify the metallic approximation. 
We consider the case of strong fields and take k ~ i/R q 
as a typical wave vector for the spatial size of the surface 
distortion. Then, using ([9]), 




2 x 10 7 



3/2 



Thus, in the fields E — 300i? c the fulfillment of inequality 
!l p > t ? _1 evidences for the favor of the metallic approx- 
imation. 



IV. THERMAL-QUANTUM CROSSOVER 
TEMPERATURE. THE DECAY RATE 

Let us turn to the thermal-quantum crossover tem- 
perature T q which separates the classical thermal activa- 
tion at T > T q from the quantum nucleation mechanism 
at lower T < T q temperatures. Here we estimate the 
thermal-quantum crossover temperature T q = T q (E) as 
a ratio of the potential barrier height to the saddle value 
S q of the Euclidean action at zero temperature. In the 
dimensional units we have then 



T q (E) = — ° 
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(in mK). 



Note that the thermal-quantum crossover temperature is 
independent of the step tension coefficient a±. This point 
is obvious since the barrier height Uq and action S q are 
both proportional to the same factor a 2 . 

At the electric fields comparable with the critical one 
E c the thermal-quantum crossover temperature, starting 
from its zero value at E — E c , proves to be extremely 
small. In the strong E 3> E c fields the thermal-quantum 
crossover temperature grows approximately as a cube of 
the field 

T q (E) w 3 x 10~ 7 (.Ey.E: c ) 3 (in mK). 

For fields E = 300i? c , we may expect a reasonable mag- 
nitude for the thermal-quantum crossover temperature of 
about 8 mK. 

Let us consider a charged crystal facet prepared in the 
metastable E > E c state with adjusting thermodynamic 
parameters such as temperature T and electric field E. 
After the lapse of some time t f, s , there will appear a 
nucleus of the rough state breaking the crystal faceting. 
Then the nucleation rate T — T(T, E) and the time of 
observation t b s are connected by the following relation 



tobs An 



1. 



where N nuc is the total number of independent nucleation 
sites and T is the nucleation rate at a single nucleation 
site. We estimate N nuc approximately as the total num- 
ber of atoms at the crystal surface, assuming that every 
atom at the surface has an equal possibility to become 
a nucleation site within the time interval t b s - For the 
crystal area of 1 cm 2 , we put 

,14 



10 1 



The nucleation rate T can approximately be estimated 



as 



r ~ j/exp(— S) 
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where v is the attempt frequency and exponent S, de- 
pending on temperature, is either Arrhenius exponent 
Uq/T or Euclidean one S q /h. The attempt frequency 
v is associated with the surface fluctuations resulting in 
nonzero bending £(t, r) of the flat crystal facet. In gen- 
eral, the frequency of crystal surface fluctuations depends 
on the magnitude of surface bending £ and the radius of 
deformation R as well. This frequency can be estimated 
by equating the kinetic energy to the potential surface 
energy in Lagrangian L fj 10|) . The order-of- magnitude 
estimate can be represented as [26] 



a>iR + ctp |C| 
PcfflCI^ 3 



1/2 



According to 27], there is one optimum path, i.e., es- 
cape path which connects the entrance point with the op- 
timum escape point and corresponds to the saddle-point 
value of the effective Euclidean action. In the quasiclas- 
sical approximation the main contribution to the decay 
rate of the metastable state is determined by such opti- 
mum escape path and its nearest vicinity. As is found 
above, at the optimum escape path a ratio of surface 
deformation £ to its radius R satisfies approximately 
|£|/-ff ~ Oii/ctQ- Then we arrive at 



of 



OoPofflCI 3 
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Next, we should estimate the equilibrium fluctuations 
of the surface bending as a function of temperature. At 
high temperatures one expects the thermal activation 
mechanism when the average energy fluctuations should 
be of the order of the temperature, i.e., aii2|C| + «oC 2 ~ 
aoC 2 ~ T. Hence, for T = 1 K, we expect 

|C| ~ 0.3 nm, R ~ 3 nm and i/~5x 10 10 Hz. 

At zero temperature the attempt frequency can be asso- 
ciated with the zero-point oscillations in the same poten- 
tial U = aiR\(\ + aoC 2 ~ aoC 2 - Using U ~ Kv{U) for an 
estimate of the ground level energy, we find 



ciy 



a h 3 p 2 cS 



1/7 



-7xl0 10 Hz, J7 = fii/~0.5K. 



Note that the magnitude of the surface bending is 
about of the interatomic spacing and the frequency has 
numerically the same order of the magnitude as the De- 
bye frequency. These magnitudes seem us reasonable. 
Thus, we have a relatively large preexponential factor 



vN n 



which can readily be compensated by macroscopically 
large potential barrier for insufficiently high density of 
charges. Eventually, if we wish to discover the process of 
the facet destruction for the time of about tens seconds, 
the exponents Uq/T or S q /h should be kept about 55. 



Due to strong exponential dependence of nucleation 
rate Y on the thermodynamic parameters T and E the 
statistical dispersion of nucleation events is not large as 
compared with the average values of the thermodynamic 
parameters at which the nucleation is mainly observed. 
The overwhelming majority of experimental points will 
concentrate in the narrow region around the average val- 
ues which correspond to the so-called rapid nucleation 
line. In essence, from the viewpoint of the time of obser- 
vation the rapid nucleation line separates the metastable 
states into two region. One region represents the long- 
living states looking as stable during the experiment and 
the other is the short-living states which decay practi- 
cally instantly. 

So, for the rapid nucleation line or the breaking field 
Eb, we may expect the following behavior. Under ther- 
mal activation mechanism at high temperatures one 
should observe 



E b {T) oc T 



' 1/4 , T>T q . 



Below the thermal-quantum crossover temperature this 
behavior should go over to the practically temperature- 
independent behavior 

E b (T) w const, T < T q . 

In the latter connection we would like to mention a 
possible effect of the energy dissipation processes. As is 
known from the quantum dynamics of first-order phase 
transitions [28| , the energy dissipation processes increase 
the effective Euclidean action and thus reduce the quan- 
tum decay rate. Accordingly, the behavior of the break- 
ing field Ef,(T) in the quantum regime should grow with 
the temperature rise and demonstrate a maximum at the 
thermal-quantum crossover temperature. However, as 
is mentioned above, the energy dissipation in superfluid 
4 He is not large at low temperatures because of negligible 
density of the normal component. That is why, we ex- 
pect only a slight manifestation of the energy dissipation 
effects in the quantum regime. 



V. SUMMARY 

To summarize, we have examined a stability of the 
charged crystal 4 He surface in the atomically smooth and 
well-faceted state below the roughening transition tem- 
perature. Like the charged crystal 4 He surface in the 
rough state, the charged crystal 4 He facet becomes un- 
stable at the same density of charges or corresponding 
critical electric field E c . However, the dynamics of the 
transition from the initial homogeneous distribution of 
charges and flat crystal surface to a spatially unhomoge- 
neous charge distribution and to a warped crystal surface 
proves to be qualitatively different. 

In the rough surface state, as the electric field ex- 
ceeds the critical value E c , the homogeneous surface state 
becomes absolutely unstable and in this sense the de- 
velopment of the charge-induced instability resembles a 



8 



second-order phase transition. In contrast, as the electric 
field exceeds the same critical value E c , the homogeneous 
state of the atomically smooth and well-faceted crystal 
surface is converted into the metastable state separated 
with a potential barrier governed by the electric field or 
charge density. The barrier height is proportional to the 
square of the linear facet step energy. 

The onset and development of the charge- induced in- 
stability at the crystal facet can be compared with the 
kinetics of first-order phase transitions accompanied by 
the nucleation and next growth of new stable phase. A 
nucleus of new phase here can be described as a fluctua- 
tion region of the crystal surface in the atomically rough 
state. The larger the charge density, the smaller the ra- 
dius of the critical nucleus. 

Unlike the charged rough crystal surface, in the electric 
fields which are tens times larger than the critical value 
E c the potential barrier still remains so high that the 



charged crystal facet, though metastable, will not break 
up for the experimentally obtainable uptime. To realize 
a breakage of the crystal facet, a few hundreds of criti- 
cal value E c should be achieved. The breaking electric 
strength Ef, depends on the temperature as well as the 
nature of the breaking mechanism. At high temperatures 
the breaking dynamics is associated with the thermal 
activation mechanism and with the quantum tunneling 
through a potential barrier at sufficiently low tempera- 
tures. 
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